Preprint typeset in JHEP style - HYPER VERSION PI-PARTPHYS-258, UMD-DOE/ER/40762-514 

Light dark matter and Z' dark force at colliders 



Haipeng An a ' 6 , Xiangdong Ji a cd , and Lian-Tao Wang 6 >f 

a Maryland Center for Fundamental Physics and Department of Physics, University of 

Maryland, College Park, Maryland 20742, USA 

b Perimeter Institute, Waterloo, Ontario N2L 2Y5, Canada 

c Shanghai Laboratory for Particle Physics and Cosmology, and Department of Physics, 
Shanghai Jiao Tong University, Shanghai 20024-0, China 

d Center for High-Energy Physics and Institute of Theoretical Physics, Peking 
University, Beijing 100871, China 

e Enrico Fermi Institute and Department of Physics, University of Chicago, Chicago, IL, 
60637 

' Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL, 60637 
E-mail: hanQperimeterinstitute . ca, xji@umd.edu, liantaow@uchicago.edu 

Abstract: Light Dark Matter, < 10 GeV, with sizable direct detection rate is an inter- 
esting and less explored scenario. Collider searches can be very powerful, such as through 
the channel in which a pair of dark matter particle are produced in association with a jet. 
It is a generic possibility that the mediator of the interaction between DM and the nucleus 
will also be accessible at the Tevatron and the LHC. Therefore, collider search of the me- 
diator can provide a more comprehensive probe of the dark matter and its interactions. In 
this article, to demonstrate the complementarity of these two approaches, we focus on the 
possibility of the mediator being a new U{\)' gauge boson, which is probably the simplest 
model which allows a large direct detection cross section for a light dark matter candidate. 
We combine searches in the monojet+MET channel and dijet resonance search for the me- 
diator. We find that for the mass of Z' between 250 GeV and 4 TeV, resonance searches at 
the colliders provide stronger constraints on this model than the monojet+MET searches. 
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1. Introduction 

The identity of the Cold Dark Matter (CDM) in the universe is one of the outstanding 
mysterious. Many possible models of dark matter have been proposed. In many promis- 
ing scenarios, such as those realizing the so-called WIMP miracle, dark matter posses 
non-vanishing couplings to the Standard Model (SM) particles. Experiments designed to 
observe dark matter through these interactions are crucial in probing the properties of the 
dark matter. 

One of such promising approaches is the direct detection, with many results being reported 
recently [1, 2, 3, 4, 5, 6]. The current direct detection of dark matter has roughly two 
frontiers. The first frontier is for dark matter mass M x in the range of 100s GeV — TeV. 
As a frequently highlighted result, the constraint on the spin-independent dark matter 
nucleon cross section, ogi, is around 10 -44 — 10 -45 cm 2 . Along this frontier, the direct 
searches are starting to set interesting constraints on the "conventional" WIMP candidate, 
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such as the LSP dark matter of supersymmetry [7, 8] or KK dark matter of models with 
extradimension(s) [9, 10]. The second frontier is in the regime of light, M x < 10 GeV, dark 
matter. In this regime, the momentum of the dark matter in the galactic halo is relative 
small, leading to suppressed recoil energy of the detector nuclei. Therefore, the bound on 
the dark matter nucleon cross section is much weaker, dropping from ctsi^IO 2 cm 2 for 
M x ~ 10 GeV to ogi^lO -39 cm 2 for M x ~ 6 GeV, for example. For the same reason, there 
is also a similar feature in the reach for spin-dependent interactions, osd- The possibility 
of light dark matter within this mass regime has received considerable attention recently 
[11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 
35, 36, 37, 38, 39, 40]. 

High energy collider experiments are crucial in the search of dark matter, independent of 
astrophysical uncertainties. Dark matter particles can be pair produced at colliders, and 
they will escape without interacting with the detector. Therefore, the basic signature of 
dark matter is jet(s) (or 7)+ fix- If dark matter is light, the production rate is deter- 
mined dominantly by the cut on $t and relatively independent of M x . Therefore, Collider 
searches can be much more powerful and provide crucial complementary information in the 
case of light dark matter where direct detection has its intrinsic limitations. 

Some model dependence in making the connection between direct detections and col- 
lider searches is unavoidable. The momentum transfer involved in the direct detection 
experiments is tiny, ^10(s) MeV. Therefore, we can typically integrate out the particle 
which mediates the interaction between DM and SM particles. The resulting operators, 
of the form Jsm • J\i provide adequate description for direct detection. The simplest ap- 
proach is to use the same operator in collider studies, treating them as contact interactions 
[41, 42, 43, 44, 45, 46, 47, 48]. At the same time, the mediator is not necessary heavy and 
therefore not in the decoupling limit. In this case, it is not a very good approximation 
to integrate out the mediator at collider energies. Such effects have already been pointed 
out and studied [41, 42, 43, 44, 45, 46, 47, 48]. Moreover, since the mediator is within the 
reach of collider energies, there are additional channels we can use to search for it, such as 
resonant production followed by the decay back to SM states. 

The main goal of this paper is to extend the earlier studies [41, 42, 43, 44, 45, 46, 47, 48] 
by including the effect and signature of a mediator which is non-decoupling. In particular, 
we combine the mono-jet + MET search with direct collider searches for the mediator, and 
demonstrate that they are both powerful and complementary. In the light of the strength 
of collider searches in the small M x and large ersi frontier, we concentrate on mediator with 
proper coupling to both SM and DM which give rise to such large cross sections. The most 
obvious example is the scenario in which dark matter is a Dirac fermion, and its interaction 
with the nucleus is spin-independent and is through a spin-1 (Z') mediator. There is a 
corresponding example with a Majorana fermion dark matter and a Z' mediator, in which 
the dominant dark matter nucleus interaction is spin dependent. 

The rest of the paper is organized as follows. In Sec. 2, we discuss the generic properties 
of Z' model with all possible renormalizable couplings to quarks and dark matter. We 
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also illustrate the transition from on-shell production of Z' to contact interactions. We 
emphasize that in order to simulate the monojet and dijet constraint correctly with a 
relatively broad resonance, one has to use kinetic width, which is discussed in detial in 
Appendix A. In Sec. 3, we discuss constraints and reaches from Tevatron and LHC in 
search for excess in the jet+ $t channel. We translate the results to constraints on direct 
detection cross sections of elastic collision between nucleon and dark matter. Results for 
both <Tsi and cjsd are presented. In Sec. 4, we discuss the dijet constraint on Z' mass and 
its couplings to the quarks. We find that if the couplings between DM and Z' is comparable 
or smaller than the couplings between quarks and Z', the constraint from dijet is stronger 
than from monojet. The constraints on couplings of the new gauge boson are shown in 
Appendix B. 

2. Leptophobic Z' model, direct detection, and searches at colliders 

If the interactions between dark matter and SM quarks are mediated by Z 1 ', including all 
renormalzable interactions, the Lagrangian can be written as 

£ = z'uKaz'q^q + gz^qi^q) + {goxi^x + go^xi^ibx)] , (2.1) 

where q and x are denoting SM quarks and dark matter particles, respectively. In general, 
Z' can also couple to leptons, and we can parameterize such extensions with B — xL type 
of couplings. In this case, the search of dilepton resonance provides the most obvious probe 
and its reach depends sensitively on x, which is not related to the direct detection of dark 
matter. In this paper, we will concentrate on the minimal case of a leptophobic Z' . 

In direct detection experiment, the momentum transfer between dark matter particle and 
the nucleus can be estimated as M^M x v cos 6 /{Mjy + M x ). Therefore, we can see that the 
momentum transfer is limited by the mass of nucleus and the speed of dark matter and 
cannot be larger than around 0.1 GeV. Therefore, in the case that Mz> is larger than a 
few GeV, the interaction between nucleus and dark matter can be described by a contact 
interaction which depends only on the Wilson coefficient gz'gD/M^,. 

We begin with a review of the connection between direct detection and various possible 
Z' interactions. After integrating out the Z' , we obtain 4 different effective operators. 
We collect them, as well as their corresponding DM-nuleon cross section and the relevant 
properties of DM- nucleus cross section, in Table 1 . It is well known that dark matter direct 
detection rate is enhanced by coherent scattering off the nucleus if the dark matter nucleon 
interaction is spin independent. Typical velocity of the dark matter in the galactic halo is 
v ~ 1CT 3 . Any process depending on the momentum exchange |Ap*Ar| = |Ap x | ~ \p x \ is 
suppressed by an additional factor of v 2 . 

In the case that Z' couples only to vector currents of both quark and dark matter fields 
(Oi), the cross section is both spin independent, and it is unsuppressed by momentum 
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Table 1: Effective operators for dark matter direct detection after integrating out the Z' , and 
the corresponding DM nucleon cross sections. We also show whether the dark matter nucleus cross 
section is independent (SI) or dependent (SD) on the nucleus spin, and whether it is dependent on 
the dark matter momentum (MD) or not (MI). Mm denotes the nucleon mass. AE is defined as 
(-^1 2~Z 9 97m 7591^0 — ^£[7^7^75?/^, where Un is the wave function of the nucleon. O3 and O4 are 
parity-odd operators inducing velocity dependent interaction between DM and nucleons. v is the 
velocity of DM in the lab frame. 



exchange. Using some typical values of gauge couplings and Mz 1 , we can estimate 

= Vf &KMl ,„y ,„y /200 GeVy 

b nM%{M N + M D ) 2 V0.3/ V0.3/ V M z> J 

Therefore, this is a very plausible way of realizing the large cross section scenario we are 
interested in probing. We are going to to mainly focus on this case in this paper. 

It is possible that 0\ is absent or strongly suppressed. An important example is when the 
dark matter particle is a Majorana fermion. In this case, dark matter can only form axial 
vector current and O2 and/or O4 describes the direct detection, with O4 typically being 
the dominant one since O2 is further suppressed by momentum exchange. We will consider 
this case separately in Sec. 3.2. 

In most of the cases considered in this paper, the U (1)' is anomalous, and spectator fermions 
need to be introduced to cancel the anomaly. In this case, the upper bound for the mass of 
spectating fermion is [49] M spectator < (Q4tt 2 / g\i)Mz> , where gz> is the coupling between 
Z 1 and SM fermions and Mz' is the mass of Z' . Potential signals from these states provide 
further information about this scenario. To focus on the most generic phenomenology, we 
will not include explicit spectators in this work. 

The size of the signal of dark matter depends on sizes of couplings go and gz> ■ In principle, 
they don't have to be related. However, if DM particle is part of the spectator, for example 
in the scenario of Ref. [50], the couplings will have to be of the same order. On the other 
hand, for example, they might be introduced through effective interactions [51] and they 
can be quite different. In this work, to illustrate the basic features of our results, we first 
assume gz> = go- From the various collider processes we get upper bounds on gz> (#d), 
with which we then set the upper bound on direct detection cross section. We then proceed 



- 4 - 



2.0 




200 400 600 800 1000 1200 1400 
M z <GeV) 

Figure 1: Cross section for monojet + MET at Tevatron as a function of Mz> with different dark 
matter masses. gz'/Mz': an d hence the direct detection cross section <7direct 5 is fixed. We assume 
9z> = go- 



to investigate how the constraints change if we relax this assumption and allow gz'/gD to 
vary. 

In collider production of dark matter particles, the relevant momentum scale, q, can be 
comparable with Mz> and the dependence of the signal rate on the model parameters are 
in general more complicated, thus the contact interaction approximation is not always a 
good one. To demonstrate this, we study the production cross section as a function of Z' 
model parameters, particularly the Z' mass, fixing the combination gz'9D/M^, = 1/A 2 
and hence the direct detection rate. We further assume gz> = go for simplicity. We have 
at parton level, 

$ZJ_ _ \ / 9 n\ 

Ccollidcr OC 2 „ — , „. 2 , {Z.6) 



{q 2 _ M | /)2 + {g2 V _ A2)2 + a2 W 



M\, " > 1 " (15hr) 

where Tz'(q 2 ) is the width of Z' with M|/ replaced by q 2 , and a is a constant depending 
on couplings and the number of degrees of freedom. 

Fig. 1 shows the monojet + MET cross section (to be discussed in detail in the next section) 
at Tevatron as a function of Mz' with A fixed to 300 GeV, for two different values of M x . 
For comparison, we have also shown the production rates predicted by using a contact 
interaction with coefficient 1/A 2 . When Mz> < 2M X , Z' cannot be on shell in this process. 
The parton level cross section for qq — > (Z'* — > xx) + X can be written as 

<5"collidcr ~ — — • (2.4) 

M^J A + a (I&F 

The cross section is suppressed by (M|,/4M 2 ) 2 at small M^. In the case that Mz> > 2M X , 
Z' can be produced on shell. Using the narrow width approximation, 

a ~ a pTod (Z') x Br(Z' -> xx) ■ (2.5) 
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Figure 2: Dominant partem level diagrams for pp — > monojet + MET. (b) is also the dominant 
parton level process for pp — > monojet + MET in LHC. 

Therefore, the threshold behavior near Mz< ~ 2M X in Fig. 1 is due to the enhancement 
from resonant production. 

For larger Mz 1 , the interaction between quarks and dark matter becomes a contact interac- 
tion which is shown by the flat tail in Fig. 1. For very large Mz>, typical q 2 is independent 
of Mz 1 and constrained by parton density to be q 2 <C Mz>- The production rate asymp- 
totes to a constant, which is less than the prediction from contact interaction due to the 
width term on the denominator in Eq. (2.3). 

In Eq. (2.3), the kinetic width of the intermediate particle is used instead of the Breit- 
Wigner approximation. In our discussions, since Z' is assumed to couple to all flavors of 
quarks universally and the gauge coupling is allowed to vary in a wide range, the width of 
Z 1 can be broad. Therefore, the Breit-Wigner approximation with a large constant width 
cannot be seen as a good approximation in this case. A detailed discussion of this effect 
can be found in Appendix A. 

3. Monojet + MET searches at the Tevatron and the LHC 

We consider dark matter particle production in association with a jet radiated from the 
initial state parton. Some examples of the Feynman diagrams are shown in Fig. 2. This 
process, arguably being the most model independent, has been studied [41, 42, 43, 44, 45, 
46, 47, 48, 52]. We present our study within our Z' framework, taking into account the 
effects of kinematic Z' width as described in Appendix A. The experimental data is from 
Tevatron and LHC searches with L ~ 1 fb -1 . The result is shown as the constraints on 
direct detection cross section, whereas the constraints on the couplings are presented in 
Appendix B. 

3.1 Constraints for the spin-independent case 

Most of the existing searches for new physics in the monojet+MET channel are carried out 
in the context of the large extra dimension (LED) model [53]. We will use those results to 
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set limits on our Z' model. Being a non-renormalizable model, the pt (and Pt) distribution 
predicted by LED is somewhat harder than the Z' model, and different cuts can be chosen 
to improve the sensitivity. 

The strategy we use to get constraint from monojet search on the Z' model follows from 
Ref. [54], by defining a x 2 > 

2 = [iVobs ~ iVsM ~ N Z >} 2 , 1 

X N z > + N SM + a 2 SM ' 1 ' J 

where -/V bs is the number of events observed in the experiment, Asm is the prediction 
from SM, cjsm is the uncertainty of the SM prediction including statistic and systematic 
uncertainties. Then we require x 2 < 2.71 to get 90% constraint on the contribution from 
the Z' model, N%i. 

The selection cuts used by various searches, which we also adopt, are listed in Table. 2. 
CDF group used 1 fb" 1 of data [55], and two sets of cut on Pt (MET) and the leading jet 
were used. The constraints on the Z' model from LowPT cut is always stronger than the 
one from HighPT cuts. Therefore, we will use the lower PT cut to set limits. With this 
set of cuts, 8449 events have been observed, which is consistent with the SM background 
prediction of 8663±332. ATLAS and CMS have published their analysis of monojet + MET 
with a luminosity of 1 fb _1 , we will concentrate on the ATLAS search in our study [56]. 
They did analysis with three different PT cuts, namely LowPT, HighPT and veryHighPT, 
as shown in Table 2. 



CDF 1 fbr 1 


p T > 80 GeV, p T (ji) >80 GeV. 
Events with additional jets are vetoed if 
Pr(h) > 30 GeV or p T (j 3 ) > 20 GeV. 


CDF 6.7 fbr 1 


p T > 60 GeV, p T (ji) >60 GeV. 

A second jet is allowed if 20 GeV < pr{h) < 30 GeV. 


ATLAS LowPT 


p T > 120 GeV, one jet p T (ji) > 120 GeV, \r}{ji)\ < 2. 
Events with a second jet are vetoed 
if Prfo) > 30 GeV and |r/(j 2 )| < 4.5. 


ATLAS HighPT 


Pt > 220 GeV, p T (ji) > 250 GeV, \ v {ji)\ < 2. 

Event is vetoed if there is a second jet with pr{j2) > 60 GeV 

or A(f>(j 2 ,$ T ) < 0.5 and \r,{j 2 )\ < 4.5. 

Any additional jet with 7/(^3) < 4.5 must have prijz) < 30 GeV. 


ATLAS vertHighPT 


p T > 300 GeV, one jet with p T {ji) > 350 GeV. 
1 77(.7i)| < 2, and events are vetoed if there is a second jet with 
|?7(72)| < 4.5 and with either prih) > 60 GeV or A(j 2 ,#r) < 0.5. 
Any further jets with 77(73) < 4.5 must have pr(j?i) < 30 GeV. 



Table 2: Different cuts defined by CDF and ATLAS in monojet + Pt searches, ji and j'2 denote 
leading and second leading jets respectively, j'3 labels any other jets. 



In our analysis, we generated events at parton level, and applies selection cuts. A full anal- 
ysis would also include further steps of parton shower and merging it with hadronization. 
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Such studies have been carried out in the case of contact operator [41, 42, 43, 44, 45, 46, 54], 
for which the results from different couplings can be obtained from simple scaling. In our 
case, we in principle have to carry out full simulations for signal at every possible mass and 
coupling in our scan . Instead, we match our parton level result in the contact-interaction 
limit (large Mz>) with the full studies in Ref. [42] (Tevatron) and Ref. [54] (LHC) and ob- 
tain a scale factor. This scale factor is then applied to our model for general Mz> and gz>- 
Although just an approximation, this does capture the main difference between a parton 
level and a full study in our case. For the CDF simulation, the cut we are using is that at 
parton level, the px of the jet should be larger than 80 GeV. For ATLAS simulation, we 
require that the px of the jet should be larger than 120 GeV, 220 GeV and 300 GeV for 
LowPT, HighPT and veryHighPT, respectively. For CDF, if the mass of DM is smaller 
than 100 GeV, the scale factor is about 0.6, about 0.4 for ATLAS LowPT cut, and about 
0.5 for HighPT and veryHighPT cuts. The difference between parton level simulation and 
the full study is mainly caused by the fact that after parton shower, the energy and px of 
the clustered just tend to be smaller than the original parton. Hence, the scale factor is 
almost always smaller than one. 
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Figure 3: Monojet constraints on direct detection cross sections for gz< = go and M x = 5 
GeV. The solid, dashed and dotted red curves arc for ATLAS Monojet constraints with VeryHighPT, 
HighPT and LowPT cuts described in Table 2. The green solid curve is the monojet constraint 
from CDF. The dashed green and blue curves are constraints from CDF and ATLAS dijet resonance 
searches. The solid blue curve is LHC 5cr reach assuming a centre-of-mass energy of 14 TeV and a 
luminosity of 100 fb _1 . 

Assuming go = 9z', we convert these constraints into limits on the direct detection cross 
section. The constraints on the spin independent cross section from Tevatron and LHC 
monojets searches [55], for fixed dark matter mass M x = 5 GeV, are shown in Figs. 3 
and 4. Very recently, this search has also been carried out using CDF Run II data set 
with a luminosity of 6.7 fb _1 [52]. A different signal region was chosen, which is also 

1 We note that there is no simple dependence of the monojet cross section on Mz' and gz' which is valid 
for the full parameter space. For example, g Z ' enters the Z' width. Therefore, the cross section is not a 
simple function oc g\, , especially when the Z' is broad. Therefore, we can not obtain correct results just by 
scaling from the simulation using a particular set of (gz 1 , Mz')- 
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Figure 4: Monojet constraints on direct detection cross sections in the case of small Mz>, assuming 
9z' = 9d and M x = 5 GeV. 

shown in Table. 2. They did a binned study in the signal region, and they translated 
their constraints on the generator level rate of the monojet + MET in the signal region for 
mediator masses of 100 GeV and 10 TeV, respectively. To incorporate it into our study 
with general mediator masses other than the two chosen values, we do a interpolation to 
get the corresponding constraints. The corresponding constraints on direct detection cross 
section is shown as the dotted black curve in Fig. 3, where we can see that the new cuts 
is different from the one set by ATLAS with VeryHighPT cuts and the previous CDF cuts 
with 1 fb _1 . 

For very large Mz>, we can effectively integrate out the Z' . The resulting contact interac- 
tion provides a good approximation even at LHC energies. In this limit, both the direct 
detection and monojet+MET cross sections depend on the same combination g^gf^/M^,, 
therefore the limits shown in Fig. 3 approach a constant value for very large Mz> ■ We can 
also see that the contact-interaction limit is reached at larger Mz' for searches at higher 
energies and more sensitive cuts, as expected. The limits become stronger for interme- 
diate values of Mz>, since in this regime, the Z' can be produced on-shell, leading to a 
significantly enhanced cross section for the monojet+MET process. When Z' mass is com- 
parable or less than the kinematical cuts used in the searches, the monojet+MET cross 
section starts to be less sensitive to M%i. In this regime, the monojet searches are effec- 
tively setting limits on g^,, while direct detection still depends on g 2 Z ig 2 D /M^,. Therefore, 
the limits becomes weaker in this range of Mz 1 , as shown in Fig. 3. The constraints for 
very small Mz> is shown in Fig. 4. We see that in this case, the constraints from collider 
searches are weak, mainly due to the M^, dependence on the direct detection cross section. 
As we will see later in this paper, the collider search to Z'-like resonances can not provide 
useful constraint in this regime either. It remains a challenge to find better probes for 
such light Z' with only hadronic decay modes. The "kink" feature in Fig. 4 is due to the 
threshold effect around the point at which 2M X > Mz>, where the signal process can only 
proceed through an off-shell Z' . 
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M x (GeV) 

Figure 5: Monojet constraints on direct detection cross section. The red, green, blue curves are 
for M z > = 100, 300, 1000 GeV, respectively, together with XENON 100 and CDMS constraints. 

Fig. 5 shows the Tevatron and LHC constraints with several fixed values of Mz>, as a 
function of dark matter mass M x . In the case of M x smaller than around 5 GeV, the 
constraints from collider physics become stronger than that from the direct detection ex- 
periments. In the regime where M x <C q, where q ~ max(M^',Py lt ) is the hardest scale in 
the process, the production rate is approximately independent of M x , and is a function of 
9z'i 9d, and Mz>- Given a set of gz', gD-, and Mz', direct detection can depend further on 
M x through the dark matter nucleon reduced mass M* = M^M x /(Mn + M x ). However, 
this dependence is rather weak for M x ~ O(10) GeV since M* ~ M^r. Taking together, we 
expect the limits derived from collider searches are rather insensitive to the dark matter 
mass My. In contrast with the steep weakening of the direct detection bound for light 
dark matter, collider searches are particularly powerful in this regime. For heavier dark 
matter, the visible "kink" -like feature around 2M X ~ Mz> in the curves are due to the 
transition from 2 — > 2 production process pp(p) — > Z'+jet followed by decay Z' — > xXi to 
2 — > 3 production process pp{p) — > (Z f * — > xx)+J e t which has a much smaller production 
rate. For example, there is such a feature on the red curve in Fig. 5 near M x ~ 60 GeV for 
Mz' = 100 GeV. One can see that this turning point also shows up at the green and blue 
dotted curves. Notice there is also a second turning point, for example ~ 300 GeV for the 
7 TeV LHC, at which the dark matter is too heavy to be produced. 

The new physics search potential in the pp — > monojet + MET channel has been studied at 
14 TeV center-of-mass energy with luminosity of 100 fb" 1 in Ref. [57]. With the requirement 
that MET is larger than 500 GeV, they found that the SM background is about B = 2 x 10 4 
events. For a 5 a discovery, we require that the signal should be larger than than 5\/~B. 
The 5 a reach is shown in Fig. 3 for a light WIMP and in Fig. 5 for Mz' fixed to be 300 and 
1000 GeV. One can see that for Mz 1 around 300 ~ 1000 GeV, LHC can potentially reach 
the interesting region where anomalies are claimed by CoGeNT [5, 58] and CRESST [6] 
direct detection. 
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Figure 6: Constraints on spin-dependent direct detection cross section between dark matter par- 
ticle and the nucleon (proton or neutron). The solid and dashed coloured curves are for ATLAS 
monojet constraint with VeryHighPT cut and CDF monojet constraint, respectively. The red, 
green, and blue curves are for Mz> = 100,300,1000 GeV respectively. The black solid curve is 
the constraint for neutron scattering from XenonlO [59], the thick dashed curves is that of proton 
scattering from SIMP2010 [60]. The thin dotted dash curves and dotted curves are from COUPP 
2010 [12] and PICASSO 2009 [61] on the constraint for proton scattering. The constraint on ctsd 
for the neutron can be got from the constraining curves scaled by a factor of (AE^/AEj,) 2 , which 
arc defined in the caption of Table 1. 

3.2 Constraint for the spin dependent case 

We have so far concentrated on the cases that Z' couples only to the vector currents of 
quarks and dark matter. If Z' couples to the axial vector current of dark matter particle, 
the direct detection cross section depends on either its momentum or the nucleon axial 
charge AS. If dark matter particle is a Majorana fermion, it can only couple through an 
axial vector current and only operators of O2 amd O4 in Table 1 are relevant for the direct 
detection. We generically expect both of these operators to be present. If Z' couples only to 
the vector current of quark fields (O2), scattering amplitude is proportional to \p\, and the 
direct detection cross section is suppressed by a factor of v 2 ~ 10~ 6 . Therefore, we focus 
on the case in which O4 dominates. The collider signal is largely insensitive to the details 
of the couplings since the typical momentum exchange there is the hard scale of scattering 
process. Assuming gz' = gD, the collider constraint on the spin-dependent is shown in 
Fig. 6. For the dark matter mass smaller than around 100 GeV, Tevatron constraint is 
already much stronger than that from direct detection. The dependence on the mass of 
mediator is similar as in the case of spin-independent scattering. The features of the curves 
are similar to the spin-independent case studied in the previous subsection. 

4. Constraint from Dijet Searches for Z'. 

If Z' is sufficiently light and its coupling to quarks is unsuppressed, direct searches for its 
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collider signal, such as a resonance in the di-jet final state, should be promising. There are 
two main approaches to search for Z' in the di-jet channel. First, one can perform a "bump 
hunting" in the dijet invariant mass distribution. One can also look for deviations in other 
jet kinematical distributions, such as angular correlation. In this section, we consider the 
collider limits on both searches and their constraint on the direct dark-matter detect cross 
section. 

4.1 Dijet resonance search with gz' = go 

CDF collaboration has studied the dijet final state constraints for Z' model with 1.13 fb _1 
data [62] . ATLAS collaboration has also published a study of dijet search of new resonances 
with a luminosity of 1 fb _1 [63]. In both studies, they fit the dijet invariant mass spectrum 
with the function f(x) = pi(l — x) P2 x p:i+P4lnx to model the QCD background, where pi are 
parameters and x = mjj/ ' \fs with rrijj the dijet invariant mass and y/s the center-of-mass 
energy of the collider. In our study, we use the Bayesian method described in Ref. [64]. 
We first assume that the prior probability density function of the coupling gz> to be a 
uniform function between and 4tt. For each value of gz r , we simulate the signal and 
calculate the likelihood function, from which the posterior probability density function is 
obtained. 

The constraints on the direct detection cross section from CDF and ATLAS, assuming 
9z> = 9d and M x = 5 GeV, are shown as the green and blue solid curves in Fig. 7. In the 
region 900 GeV < Mz' < 4 TeV, dijet resonance search in ATLAS gives the most stringent 
constraint, while for Mz> around 300 GeV to 900 GeV, the CDF dijet resonance search 
gives the most stringent constraint. 




M T (GeV) 

Figure 7: Constraints on spin-independent direct detection cross sections assuming gz* = 9d 
and M x = 5 GeV. The solid red curve is the constraint from ATLAS monojet constraints with 
VeryHighPT cuts described in Table 2. The green solid curve is the bound from di-jet resonance 
search by CDF. The dashed green and red curves are constraints from CMS and DO dijet angular 
distribution measurements. The solid blue curve is the constraint from ATLAS di-jet resonance 
search with 1 fb _1 data set. 
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The corresponding constraints in the case of heavier dark matter can also be obtained by 
simple scaling of the constraints presented in Fig. 7. In the case of go = gz', the constraint 
from dijet resonance searches depends very weakly on the mass of DM since the Z' — > XX 
decay channel only contributes ~ 5% of the width of Z' . Therefore, the threshold effect on 
Tz' at Mz> ~ 2M X is negligible, and the constraints on Z' mass and couplings from CDF 
and ATLAS dijet resonant searches are also applicable to the case of heavy DM. Therefore, 
for M x ^ 5 GeV, one could just scale the di-jet constraints from Fig. 7 by a factor of 

M 2 (M N +5) 2 
5 2 (M N +M X )' 2 • 

As shown in Fig. 7, In the region between 200 GeV < Mz> < 4 TeV, the dijet resonance 
searches provide significantly stronger constraints than the mono-jet + MET searches. 

4.2 g z > / go 

In some sense, the assumption gn = gz> is overly simplistic in the dijet case, because at 
large Mz 1 , the dijet cross section constraints only gz 1 with little effect on g^. In fact, 
the dijet constraint is so strong that go can have a very large value before runs into 
contraction with monojet constraint. A large go certainly weakens the direct detection 
constraint considerably. 

Let us study the effect of allowing gz> ^ gr>- The constraints for go/gz' = 1,3,5,10,20 
for monojet and dijet resonance studies are shown in Fig. 8. The dependence of dijet cross 
section on gp is only through next leading order effects, such as the Z' width. Therefore, 
we expect the dijet constraint on gz 1 depends on gp weakly unless the invisible width of Z' 
dominates over the visible in which case the constraint on gz> becomes weaker. Therefore, 
we can see that the constraint on direct detection cross section still scales roughly as 
(or slightly stronger), and it becomes much weaker in the case of go 3> gz 1 - On the other 
hand, the monojet cross section is proportional to gz'9o unless decay channel Z' — > XX 
dominates. Therefore, one expects that the constraint on direct detection cross section 
does not change much with gp. From Fig. 8, we see that in the region go i?> gz', monojet 
constraints may dominate over dijet constraints. 

4.3 Dijet Angular Distribution 

CDF, DO, ATLAS and CMS collaborations have done searches for quark compositeness 
through study of dijet angular distribution [65, 66, 67, 68]. The quantities of interest is 
the normalized differential cross section a = (l/CTdijet)(dadijet/dXdijet), where Xdijet = (1 + 
| cos (9*|)/(1 — | cos 8*\) and 8* is a scattering angle in the rest frame of the two partons. 

In our study, the QCD background is simulated using Pythia 8.1.4.5 [69], and Fastjet [70]. 
We use AntiKT jet algorithm with R = 0.4. The new physics contribution is simulated 
using a private code [71]. For each group of dijet centre-of-mass energy defined in Refs. [66, 
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Figure 8: Comparing monojet and dijet constraints. The solid, dashed and dotted curves are for 
ATLAS dijet resonance search, ATLAS monojet search with VeryHighPT cut and CDF dijet search, 
respectively. The red, green, blue, pink and black are for go/gz' — L 3, 5, 10, 20, respectively. The 
mass of DM is assumed to be 5 GeV. 



68], we calculate the x 2 which defined as 

/=ncw , -QCD = exp\2 
j °exp ' °QCD 

where af ew , cr^ CD and cr^ xp are the new contributions, QCD background and experimental 
value in the i-th bin for certain Mjj group, respectively. <5 CX p and <5qcd are the uncertainties 
of experimental values and QCD background. To get 95% C.L. constraint on gz> for certain 
values of go and M^, we require that in each rrijj group the possibility to get calculated x 2 
should be smaller than 0.05. The constraints on gz> from CMS and DO are shown in Fig. 9, 
where the red and green curves are for DO and CMS respectively; and the corresponding 
constraints on direct detection cross sections are shown in Fig. 7. 

Since Tevatron is a pp collider, the main background is from qq — > jj and gg — > jj. The 
dominant contribution to the signal is from qq — > Z' — > qq, where Z' can be either on or 
off shell, gg — > gg provides dominant background in the energy region of V~§ < 300 GeV. 
However, it drops steeply at vl ~ 500 GeV, where qq — > jj becomes dominant with a much 
smaller rate. At the same time, Z' with Mz' ~ 500 GeV can still be produced on-shell. 
Therefore, we see from red curve in Fig. 9 that the constraint gets stronger at around 500 
to 800 GeV. For larger Mz>, Z' on-shell production is strongly suppressed by the steeply 
falling PDF. As a result, the constraint on the coupling gets weaker and eventually reaches 
the limit of the contact interaction, which is illustrated by the plateau of the red dashed 
curve in Fig. 7. The height of the plateau can be interpreted as A ~ 2 TeV for a quark 
composite operator (27r/A 2 )(g7^<7) 2 which agrees with the result from the compositeness 
search at DO [66]. 

At the LHC, the major background comes from gg — > jj and qq — > jj. The signal contains 
two contributions which are shown in Fig. 10, where (a) is an 1/Nq suppressed interference 
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Figure 9: The constraints on the Z' coupling, gz>, from di-jet angular distribution measurements 
at DO (red) and CMS (green). 




Figure 10: Dominant contributions to dijet angular distribution at LHC. 



between QCD and the new physics which is enhanced by the parton distribution function of 
the valence quarks; and (b) is suppressed by parton distribution function of q but will still 
be significant if Z' can be produced on shell. It turns out that in the case of light Z' with a 
mass smaller than around 500 GeV, the contribution from qq initial state dominates. Since 
the first energy bin for this study in CMS starts from m^ in = 250 GeV, therefore we expect 
the resonance enhanced qq contribution to the angular distribution can only be significant 
when Mz 1 is comparable with 7n^ in . This is the reason for the local enhancement of the 
constraint around Mz> ~ 250 GeV, shown in the green curve in Fig. 9. There is no such 
a dent in the red curve since for the DO study, the constraint is always from the last few 
bins of the dijet center of mass energy. 

The constraint on gz' from dijet angular distributions is not as sensitive to the width of 
the mediator as in the case of dijet resonance search. Therefore, the constraint on direct 
detection cross section can be estimated by simply scaling the curves in Fig. 9 by a factor 

of 9d/9z" 

From Fig. 7, we can see that the current constraints from dijet angular distribution search 
is not as strong as the one from dijet resonant search. In longer term, we expect angular 
correlation measurement will be more relevant for the case of very large Mz' ■ 



5. Conclusion 



In this work, we consider the possibility of probing light dark matter at the colliders. We 
interpret the current limits and potential search reach at the colliders as constraints on 
the direct detection cross section. As already been demonstrated in earlier works, such an 
approach can yield stronger, albeit model dependent, constraint in comparison with the 
direct detection for the light dark matter M x ~ 5 GeV. The simplest approach to connect 
the collider and direct detection limits is to work in the limit in which the mediator is so 
heavy that it can be integrated out both for direct detection and for collider searches at the 
Tevatron and the LHC. Focusing the monojet+MET observable, the connection between 
these two set of observables can be established in a straightforward manner. 

We argue that while using contact operators is probably the most straightforward approach 
with the most direct connection between dark matter direct detection and collider searches, 
it is important to consider the case in which the mediator of dark matter and SM interac- 
tions will also be accessible at high energy colliders. Earlier studies have emphasized the 
effect of light mediator on the monojet + MET channel. In this work, we adopt the strat- 
egy of combining the search in direct dark matter production channels of mono-jet (photon) 
+ MET, and the search of mediator directly at the colliders. We demonstrate that the two 
approaches are complementary. 

As a concrete example, we study the case in which the dark matter is a Dirac fermion with 
its interaction with the SM mediated by a spin-1 Z' through vector like couplings. Such a 
Z' can be searched directly in dijet channels. We studied the constraints on the Z' model 
from both dijet and mono-jet + MET searches at the Tevatron and the LHC. In the case of 
a strongly coupled Z\ we have carefully taken into the effect of the finite resonance width. 
We then combined the constraints in dijet resonance searches with the mono-jet + MET 
searches, and interpreted them as constraints on the direct detection cross sections. In the 
case of spin-independent interaction, the monojet constraint is stronger than the direct 
detection constraints only in the case that M x is smaller than about 5 GeV. At the same 
time, if the mass of the mediator is larger than around 250 GeV, the constraint from dijet 
search is stronger than the monojet constraints. In the case of spin-dependent interaction, 
the collider constraints become much stronger then the direct detection constraints. 

The process of mono-photon + MET is also studied in CDF [55] and a theoretical study has 
been done for LHC in Ref . [54] , which shows that the mono-photon constraints are relatively 
weaker than the mono-jet constraints due to that the rate is much smaller although the 
mono-photon signal is cleaner than monojet. 
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A. Propagator with large width 

In the prototype model described in this work, Z' is assumed to couple all the quarks 
universally. Therefore, if gz> becomes large, Z' becomes a broad resonance and the Breight- 
Wigner approximation can no longer be used as a good estimation. In the broad resonance 
case, the kinetic dependence should be included in the propagator, which can be written 

as 

2. * AT M ' 

s -m z + i^sT tot {s) 

where i\/sTt t(s) includes all imaginary parts in the two-point correlation function of Z' . 
The diagrams are shown in Fig. 11. The diagrams can be separated into two groups, in 
the first group shown in the first line of Fig. 11, it is easy to see that the each diagram is 
proportional to (g 2 N) r , where g is the coupling between the intermediate particle and Z' , 
N is the number of fermion degrees of freedom and r is the number of fermion loops. In 
our case, since the Z' universally couples to all quarks, there are at least 15 light degrees of 
fermions. Therefore, one can see that the width of Z' becomes broad when g is order one. 
And the summation of the first line in Fig. 11 dominates over others. After summing up 
all the leading order diagrams, the propagator can be written as in Eq. (A.l). In the case 
that yfs S> mf where mf is the mass of the fermion in the loop, Tz'(s) is proportional to 
y/s and can be estimated as 

IWW - j^T z , , (A.2) 
where Yz> is the physical width of Z' . 

Fixed width is used in most of packages of simulation which might not be very proper in 
dealing with broad resonances. A comparison between propagators with kinetic width and 
fixed width are shown in Fig. 12, where the solid and dashed curves are for the kinetic and 
fixed widths cases, respectively. Different colors are for different choices of the 't Hooft 
coupling defined as g\/N changing from 2 to 10. One can see that in the case of g\/~N = 10, 
the deviation of the fixed width one from the kinetic one can be as large as a factor of ten. 
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In our simulation, we modified the CalcHEP 2.5.7 under the instruction the authors 2 so 
that the effect of kinetic width can be realized. 




Figure 11: Diagrams of Z' two-point correlation function. 
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Figure 12: Comparison between the kinetic propagator (solid curves) and the fixed width ones 
(dashed). Red, green, blue, purple and black curves are for g* — 2,4,6,8,10, respectively. The 

curves are showing the absolute square of the propagator which is ^ — ^ • 
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B. Constraints on couplings of Z' 

We present the constraints on gz> from monojet and dijet experiments. In Fig. 13, we show 
the red and green curves corresponding to monojet-only constraints from CDF and ATLAS 
with VeryHighPT cut, respectively, asssuming M x = 5 GeV and gz> = 9d- We see that for 
M z i smaller than 2M X , the constraint is around 0.4. For Mz' between 2M X and 100 GeV, 
because of the resonance production, the constraint is much stronger. For Mz> > 200 
GeV, the constraint is in the effective theory form gz'9D/Mg,, or gz 1 is linear in Mz 1 - 
Apparently, ATLAS constraint is better in the high-M^/ region. For the same low M x , the 
dependence on the ratio of gol 9z' is shown in Fig. 14 where the solid and dotted curves are 
now constraints from dijet studies of ATLAS and CDF, whereas the dashed curves show 
the constraints from ATLAS monojet study with VeryHighPT cut. The red, green, blue, 
purple and black colors are corresponding to goldz 1 = 1,3,5,10,20, respectively. Two 
important features can be commented on. First, at large Mz>, the dijet contraint is much 

2 We thank A. Pukhov for help in modifying CalcHEP. In particular, the relevant modifications are made 
in the file sqme_aux.inc in the folder c_source/sqme_aux to invoke the kinetic width. 
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Figure 13: Constraints on gz> from monojet studies at CDF (red curve) and ATLAS (green curve), 
assuming M x = 5 GeV and gz> = <3d- 
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Figure 14: Monojet and dijet constraints on gz' for gojgz' — 1,3,5,10,20, which are labeled 
by red, green, blue, purple, black curves respectively. The solid, dotted and dashed curves are 
for ATLAS dijet resonance search, CDF dijet resonance search and ATLAS monojet search with 
VeryHighPT cut, respectively. We fix M x = 5 GeV. 

stronger than monojet, particularly from LHC. At very large Mz>, the dijet constraint is 
mostly on gz' independent of go- As Mz 1 becomes smaller, gp affects the width of Mz>, 
larger gp weakens the constraint on gz>- However, to weaken this constraint to the same 
level as the monojet, one would allow an extremely large gn, beyond the effective range 
of the Z' model. Thus the tightest constraint happens for gz 1 from dijet and gu from a 
reasonable theoretical bound. Second, at Mz> less than few hundred GeV, the constraint 
is mostly from monojet, and the relevant parameter is just the product godz 1 ■ All results 
are nearly independent of dark matter mass so long as M x < 50 GeV or so. 
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